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A quark-model framework for studying nuclear medium effects on nucleon resonances is described 
and applied here to pion photoproduction on the deuteron, which is the simplest composite nucleon 
system and serves as a first test case. Pion photoproduction on nuclei is discussed within a chiral 
constituent quark model in which the quark degrees of freedom are explicitly introduced through an 
effective chiral Lagrangian for the quark-pseudoscalar-meson coupling. The advantage of this model 
is that a complete set of nucleon resonances can be systematically included with a limited number of 
parameters. Also, the systematic description of the nucleon and its resonances at quark level allows 
us to self-consistently relate the nuclear medium's influence on the baryon properties to the intrinsic 
dynamic aspects of the baryons. As the simplest composite nucleus, the deuteron represents the first 
application of this effective theory for meson photoproduction on light nuclei. The influence of the 
medium on the transition operators for a free nucleon is investigated in the Delta resonance region. 
No evidence is found for a change of the Delta properties in the pion photoproduction reaction on 
the deuteron since the nuclear medium here involves just one other nucleon and the low binding 
energy implies low nuclear density. However, we show that the reaction mechanism is in principle 
6JT), sensitive to changes of Delta properties that would be produced by the denser nuclear medium of 

^ ■ heavier nuclei through the modification of the quark model parameters. 
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I. INTRODUCTION 

The question of the influence of the nuclear medium on the properties of the nucleon and its resonances has attracted 
a lot of attention recently both in experiment and theory. It has also been one of the challenges to understanding 
strong QCD 0,0 where the connection between the nucleon internal degrees of freedom and nuclear medium effects 
has been a key issue. Historically, nucleon resonances have been among the richest sources of information on low 
energy QCD phenomena. In particular, nucleon resonance excitations by electromagnetic probes and their decays via 
the strong interaction have been an ideal method for such a study. In experiment, the advent of high intensity photon 
beams at JLab, ESRF, ELSA, and MAMI has greatly improved the data for pion production on single nucleons at 
the resonance regions over the last decade. These precise measurements significantly improved our know ledg e about 
nucleon resonances and further constrained theoretical phenomenologies EJIM II S U, 13 El 13 13 111 13 El 
However, our knowledge of resonance properties in nuclear reactions [^TilTl^ is still very limited, although great 
strides have been made by several groups to tackle this issue with various prescriptions |23l 123. |25| . Because of the 
complexity of the strong interaction in nuclear reactions, many phenomena may contribute, making the theoretical 
analyses difficult. For instance, recent experimental work [2(j on 7 + d — > n° + X and 7 + Ca — > tt° + X suggested 



^ \ a controversial role for the Z?i3(1520) by an empirical property change of the resonance |25l l27j. It is essential for 
any theoretical model to introduce as few parameters as possible. This is our motivation for providing a quark-level 
prescription as a first step in this direction. 

We propose a quark model approach for meson photoproduction on light nuclei, based on the recent success of a 
chiral effective theory for the pseudoscalar meson (jr, rj, K) photo- and electroproduction |l5lll6lll7lll8lll9j . Starting 
with an effective chiral Lagrangian for quark-pseudoscalar-meson coupling, this framework would in principle permit 
a self-consistent inclusion of all the s- and w-channel resonances. This is a sensible approach to the investigation of 
the effects of the nuclear medium on the resonances since we believe that such effects should have an overall influence 
on all the active baryons, i.e. both the nucleon and its resonances. 

In this work, we will re-visit — > n°d. The deuteron, as the simplest composite nucleon system, provides the 
simplest example for which the effective theory can be directly applied. We will then explore the response of baryons 
(described by the quark model) to the nuclear medium, and establish the relation between nuclear medium effects and 
baryon intrinsic dynamics. Of course while we acknowledge that medium effects are hardly likely to have a significant 
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effect for the case of photoproduction on the deuteron, the simplicity of the system in this preliminary work allows 
for a systematic study of the way such medium effects will be manifest in heavier systems, the subject of future work. 

It will be helpful to briefly review the theoretical approaches to — > ir°d available in the literature. The wi dely 
applied isobaric model has been successful in giving a reasonable description of the pion cross sections (see e.g. Ref. 1281 
[23)- Phenomenological methods have also been carried out based on analyses of single nucleon interactions |30L l3l( . 
A relativistic quantum hadronic dynamical approach has also produced similar results (s^. More recently, a coupled 
channel approach adopting the Blomqvist-Laget operators was developed [3l| • In general, these approaches have 
reproduced the experimental data reasonably well in the Delta resonance regions due to the dominance of the Delta 
magnetic dipole. 

A key question in nuclear reactions is the interplay of ambiguities from the reaction mechanism and short-distance 
components of the nuclear wavefunctions. Promisingly, for the two nucleon system, such a problem can be clarified due 
to the sensitivities of transition mechanisms to certain kinematics. This feature simplifies the problem and highlights 
the physics of interest. Here, we summarize several points which will support the simplifications, and enable us 
to derive novel information about the nucleon and its resonances in — > n d. Firstly, concerning the reaction 
mechanism, we neglect the rescattering contributions at the Delta resonance region due to the dominance of the 
direct scattering process above the production threshold. This is consistent with the findings of other models |29ll33|. 
Secondly, concerning the nuclear wavefunction, we are interested in the role played by the short-distance components 
and their effects on the cross section. We will show that, in the A resonance region, the deuteron £>-statc generally 
gives corrections to the backward angle cross sections, while the 5*-state dominates at forward angles. We will then 
investigate the impact of the nuclear medium effects on both of the kinematical regions. 



II. THE MODEL FOR PION PHOTOPRODUCTION 



Let us start with an effective Lagrangian for the quark-Goldstone-boson interaction yj] 

£*ff = + V +75^) ~ m]i> + ■■; (1) 

where V and denote the vector and axial currents which have the following expressions: 

K = (2) 



and the chiral transformation is, 



€ = e *»//» (3) 



where f m is the decay constant of the meson </> m , which is the Goldstone boson field. In the SU(3) symmetry limit, 
the quark field tp is 

ip(u) 

4, = I m 1 , (4) 



and the meson field 4> m is a 3(g>3 matrix: 
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(5) 



Thus, the Lagrangian in Eq. Q is invariant under the chiral transformation. 

Expanding the field £ in Eq. © in terms of <j> m , i.e. £ = 1 + i<f> m /fm + • ■ "j we obtain the standard quark-meson 
pseudovector coupling at tree level [TfiLUfij: 

H m = -^-ip-f^d^^mip , (6) 
Jm 

where the Goldstone boson field <\> m is now presented by the SU(3) octet mesons (71", 77 and K). 
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The quark-photon electromagnetic coupling is 

Hem = -eg^A^ik, r)V>, (7) 

where the photon has three momentum k, and the constituent quark carries a charge e q . 

This approach naturally places a constraint on the relative phases of those octet mesons in the SU(3) flavor 
symmetry limit, and relates the quark operators to the hadronic ones via explicit quark model wavefunctions for the 
baryons states. In this sense, the quark model transition amplitudes can be compared with hadronic ones, from which 
information about the baryon structure can be learned. 

A systematic investigation of the four channels of pion photoproduction on a bare nucleon was carried out in |16| , 
where it was shown that, up to 500 MeV photon energies, all the available observables could be reasonably accounted 
for with only one parameter for the Delta resonance coupling. Great effort was devoted in that work to restrict the 
number of parameters by taking advantage of the quark model symmetry in order that a complete set of resonances 
could be systematically taken into account. Overall quantitative agreement of the model with experimental data can 
be seen in |l6j |. In the present study, however, due to the fast fall-off of the deuteron momentum space wavefunction 
at large momenta and the photon energies considered (F 7 ~ 300 MeV), contributions from resonances in the second 
resonance regions, (Pn(1440), ,511(1535), £>i 3 (1520), 531(1620), L> 33 (1700), P 33 (1600), Pi 3 (1720), Fi 5 (1680) [H, 
etc) are generally negligible. In the first resonance region, where the Delta and Born term are dominant, it is a good 
approximation to include only these two contributions in the direct scattering process. 



III. PHOTOPRODUCTION ON THE DEUTERON 

In the 7-d cm. system the unpolarized cross section for jd — ► ir°d can be written as 

=64^2L7Tl £ \M fi (m h m f ,m, = +l)\>, (8) 

where k is the photon momentum defining the z-axis, and s is the squared total cm. energy; = 1 is the deuteron 
total spin, while m„ m/, and m 7 denote the spin projections for the initial, final state deuterons, and the incoming 
photon, respectively. The amplitudes with m 7 = —1 and +1 are not independent due to parity conservation. 
The invariant amplitude Ai fi in the direct scattering (Fig. can be expressed as 

M fi = J dpi ^(pi+q/2) (p' 2 ,q|f 7jv ^ 0A ,|p 2 ,k,) ^(pi+k/2), (9) 

where denotes the deuteron ground state wavefunction with spin projection m; T^ N _,^o N is the transition operator 
for the single nucleon interaction derived within the quark model, and is the nucleon photoproduction operator 
averaged over the isospin wavefunction of the deuteron. The kincmatical variables are defined as, P2 = — k — pi, and 

p 2 = q Pi- 

The deuteron wavefunction is defined as 

<(P)= E i L UL(\p\)Y LA (p)(LA,lm s \lm) X d lms7 (10) 

L=0,2;A,ro s 

where L = and 2 denote the S and D-state components and the spin wavefunction is 

Xim a = 'Yj (\m 1 ^m 2 \lm s )xi mi XL m2 , (11) 

rni .mi 

where x is the nucleon spin wavefunction, and mi and mi are spin projections for the spectator and active nucleon. 

The well-established Paris model ;36j wavefunction for the deuteron is adopted, which means ambiguities from the 
nuclear wavefunctions are minimised in this examination. The Paris potential is obtained from a phenomenological 
momentum-dependent meson exchange model for the nucleon-nucleon (N-N) interaction with parameters chosen to 
fit n-p scattering data for energies from the deuteron binding energy up to several hundred MeV, thus giving a reliable 
description of the deuteron wavefunction down to short n-p separations. 

To check the sensitivity of the reaction mechanism to the nuclear structure, we also carry out the calculation 
with a Hulthen wavefunction, for which the detailed description and parameters are given in Ref . 37] . This simpler 
wavefunction is generated by the Yamaguchi-Yamaguchi separable n-p interaction |38j | which has the useful feature 
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that the S and D components of the wavefunction have simple analytical forms in momentum space. The n-p potential 
in this model has parameters adjusted to give the experimental deuteron binding energy, quadrupole moment, and 
n-p triplet state effective range, and scattering length, all to high accuracy. However, it is not expected to be as 
accurate as the Paris wavefunction at larger momenta. Thus, since the long-distance asymptotic behavior of the 
deuteron wavefunction is well-known, the contrast will further pin down the nuclear structure effects arising from 
short-distance components. We neglect the D-state to D-state transition term in the amplitude in Eq.(4) but retain 
the S to D transition. The D-D transition may become important at high momentum transfers and show up in the 
backward-angle cross sections. However, this is not the focus in this work. 

Firstly, we examine the single nucleon reaction operators in — > n°d by applying the transition operators derived in 
the free nucleon reaction to the direct scattering process (Fig.[lJ. Note that the quark model parameters are determined 
in the free nucleon reactions while the deuteron wavefunction used is from independent studies. Our calculations are 
therefore parameter-free predictions of the effective theory in the quark model framework. It should also be pointed 
out that since the final state pions are treated as point-like particles, the size effects and possible medium effects on the 
mesons are absorbed into the quark model parameters. This is different from isobaric prescriptions, where the meson 
and baryon degrees of freedom could be influenced by the nuclear medium separately (see e.g. Refs. [23, 13 In 
this study, we will present a new dynamic view of this phenomenon by relating the nuclear medium's influence on the 
baryon properties to the intrinsic dynamic aspects of baryons via the quark model parameters. We note that similar 
ideas are also proposed in Ref. |39j . 

In Fig. [3 the results for the differential cross sections at E~ = 300 MeV with both wavefunctions are presented. As 
shown by the solid (Paris wavefunction) and dashed curve (Hulthcn wavefunction) in Fig. |2 these two wavefunctions 
produce the same cross sections at forward angles due to the same asymptotic behavior, while discrepancies arise at 
large angles due to differences between these two models at short distances. Although the experimental data |4£j are 
still sparse, the cross sections using both wavefunctions are in good agreement with the data over a wide angle region. 

We also present the results without the D-state contributions (the dotted and dot-dashed curve stand for the 5-state 
Paris and Hulthen wavefunctions, respectively), where the large- angle cross section differences highlight the D-state 
contributions. 

The above result is a direct examination of the quark model treatment for the elementary process, which deserves 
more careful consideration. Since nucleons in the deuteron are weakly bound, one would expect a weak effect from the 
nuclear medium. The above result indeed suggests that the baryons behave as they do in the single nucleon reactions. 
Therefore, in principle no parameter ch ang e to the baryon properties is needed here. This conclusion is consistent 
with other previous studies [2^, |2^, |3(], EI |3^, HI] . However, where our approach has an advantage is that it also 
allows us to examine the impact of the nuclear medium more systematically. We consider a mechanism that changes 
the properties of the Delta resonance due to the change of quark model parameters, e.g. quark potential strength 
within the baryons. 



IV. MECHANISM FOR MEDIUM MODIFICATION 



The quark model potential strength ah and constituent quark mass m q are explicitly related to the mass of the 
baryons in the nonrelativistic constituent quark model (NRCQM), where the nucleon and Delta resonance belong to 
the same representation |4lj| . A simple mass relation from the quark hyperfine interaction will allow us to make an 
empirical connection: 

M N (a h ,m q ) = M (a h , m q ) - ^S q (a h , m q ) , 



M A (a hl m q ) = M (a h ,m q ) + -6 q (a h ,m q ) , (12) 



where M is the degenerate mass between the nucleon and Delta before the hyperfine splitting, and Mjy and are 
the physical masses of the nucleon and Delta. The quantity 5 q is determined by the contact coupling: 

«. - £A . («) 

Ay 2irm~ 

where a s is the electromagnetic coupling. It is reasonable to assume, for the case of the deuteron, that the nucleon 
mass M/v is fixed for any set of parameters. We can then express the mass formula as: 

MA(a h ,m q ) = (^j [M A {a h ,m q ) - M N ] + M N , (14) 
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where the "barred" quantities refer to the values given by the medium-modified parameters. 

The Delta's strong decay width for A — > ir°N can be also related to the quark model parameters: 



r A 



91na\ h\ 3 (E f + M 



4tt J 8M A M^ \E f + M N 
x e - q2/3Q * , (15) 



where Ef = (M 2 ^ + Iql 2 ) 1 / 2 and co m — (m 2 + |q| 2 ) 1//2 are the nucleon and pion energies in the Delta cm. system. 
These variables are determined by the Delta mass, which hence are related to Eq. (|14|) . Also, the change of the 
potential strength will influence the form factors for baryons, i.e. both nucleon and Delta. Therefore, even though we 
fix the mass of the nucleon, its size may change due to the change of quark potential strength in the nuclear medium. 
For the Delta, we can simply express r A = r A (a/j, m q ). The irNA coupling g^NA = 9-kNnC-kN & is determined in 
the single free nucleon reaction |l6j, where C^na is the strength given by the experimental data. 

The simple relations from the static properties of baryons introduces new aspects in the reaction mechanism. If 
the medium influence leads to a significant change to and m q , from Eq. I|14f> and 1|15|) . the mass shift and the 
width change of the Delta will be correlated. As a sensitivity test, we assume the quark mass is fixed, and find that 
a change of the quark potential strength from ah — 330.31 MeV (determined in the free nucleon reaction 0|) to 
a/j = 342 MeV results in a mass increase of about 36 MeV as suggested by other calculations and total width change 
from T A = 120 MeV to f A = 156.7 MeV. 

In Fig. we present the calculated results (dashed curve) for the differential cross sections with such a property 
change to the baryons. By comparison with the unchanged-parameter calculation (solid curve) , we see that the cross 
section is very sensitive to the quark model parameters at forward angles. This can be understood as follows. The 
overall shape of the cross section is governed by the behaviour of the deuteron wavefunctions, which ensure that the 
integrand in Eq. @ is only non-negligible when the arguments of both wavefunctions are small. This is when pi is 
of the same order as both k/2 and q/2, and when k and q are parallel (forward scattering). In the cm. frame we 
also have Pi + P2 = — k and hence at forward angles where pi is large, P2, the momentum of the active nucleon, is at 
its smallest; at larger scattering angles it carries away more momentum from its interaction with the photon. Since 
the photon energy is chosen to be close to the Delta threshold in the photon-active nucleon lab system, the forward 
angles correspond to the kinematics where the Delta will contribute significantly to the cross section, and medium 
effects on the Delta will also show up most here. At backward angles, most of the momenta carried by the deuteron 
will be on the active nucleon, resulting in a shift of the photon-active-nucleon cm. energy far away from the Delta 
mass position leading to the small effects from the parameter change. 

The fact that the overall shape of the cross section is governed by the behaviour of the deuteron wavefunctions, also 
accounts for the negligible effects from the energy conservation breaking at the initial photon-nucleon and final meson- 
nucleon interaction vertices in this reaction. For denser nuclear medium, such effects may become important due to 
larger contributions from short-range components of nuclear wavefunctions. This could be a possible mechanism by 
which the nuclear medium can change the quark model parameters. 

In Fig. 01 the medium influence on the Delta magnetic dipole is compared with the free nucleon case in the 
photon-nucleon cm. frame. The largest effects are from the region where the Delta is excited close to its mass 
shell. As discussed above, for a photon energy E 1 — 300 MeV, the close-to-threshold region corresponds to the 
forward scatterings, with the off-shell region to the backward angles. When the Delta goes far off-shell, the parameter 
change effects become negligible. Thus, one would expect that the forward angle scattering is the ideal place for the 
investigation of the Delta property changes due to the medium. 

It is noted that the major contributions to the small angle cross section come from deuteron configurations in which 
the neutron-proton (n-p) separation is larger than the range of strong interactions. For these configurations any 
medium modifications to ah will be negligible. Figure [21 should therefore only be used as an indication of the relative 
sensitivity of the pion production process at different momentum transfers to a 'fixed' density-independent change 
to ah- Any density dependent mechanism which produced a change at small n-p separations (and large momentum 
transfer) of the magnitude used in our sensitivity test would be expected to be associated with much smaller changes 
at large n-p separations and small momentum transfers. 



V. SUMMARY 



We have studied the photoproduction of ir on the deuteron at the energy of the Delta excitation in a quark model 
framework with an effective chiral Lagrangian. The advantage of this approach has been that the nucleon and its 
resonances were introduced into the nuclear reaction with only a limited number of parameters [l6j . In this approach 
we propose that the nuclear medium effects on the baryon properties can be recognized through the change of quark 
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model parameters for baryons. In particular, the nuclear medium may influence the quark potential strengths of a 
baryon state, and further lead to changes to its mass and width. It thus naturally relates the resonance property 
changes to their intrinsic dynamical aspects. Mechanisms by which the nuclear medium changes the quark model 
parameters need to be addressed by investigating reactions on heavier and denser nuclei. 

We expect that a systematic investigation of reactions on heavier nuclei will establish the medium-density- 
dependence of the quark model parameters. It is also important to extend this approach to the second resonance 
region, where higher resonances such as Pn(1440), I?i3(1520), <Sii(1535), etc., can be investigated. Consequent work 
will be reported elsewhere. 
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FIG. 1: Schematic diagrams for: (a) the direct scattering process in jd — > n°d; (b) the s-channel process for the pion production; 
and (c) the w-channel process for the pion production. In (b) and (c), the photon and pion will couple to the constituent quarks 
of the baryons as illustrated in Ref. W]. 
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FIG. 2: Angular distributions for -yd — » Tv°d at E- t — 300 MeV. The solid and dashed curve denotes the calculations adopting 
the Paris and Hulthen wavefunctions with both S and D states, while the dotted and dot-dashed curves denote those with only 
the S state components, respectively. Data are from Ref. |4fl|. 
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FIG. 3: Estimation of effects from the Delta property change. The solid curve is the result of the unchanged-parameter 
calculation (the same as the solid one in Fig. The dashed curve denotes the results for a changed Delta with a flattened 
width and larger mass. Data are the same as in Fig. [21 
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FIG. 4: The Delta magnetic dipole plotted in term of the photon energy in the jN lab system. Data are from SAID 
analyses (4^. 



